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THE DISK AROUND COKU TAU/4: CIRCUMBINARY, NOT TRANSITIONAL 1 
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ABSTRACT 

CoKu Tau/4 has been labeled as one of the very few known transition disk objects: disks around 
young stars that have their inner disks cleared of dust, arguably due to planetary formation. We 
report aperture-masking interferometry and adaptive optics imaging observations showing that CoKu 
Tau/4 is in fact a near-equal binary star of projected separation ^53 mas (~8AU). The spectral 
energy distribution of the disk is then naturally explained by inner truncation of the disk through 
gravitational interactions with the binary star system. We discuss the possibility that such "unseen" 
binary companions could cause other circumbinary disks to be labeled as transitional. 
Subject headings: stars: low-mass, brown dwarfs 
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1. INTRODUCTION 

The lifetime of disks around young stars sets the time- 
scale for giant planet formation. Most young stars have 
disks that are optically thick at all wavelengths or have 
no observable disk at all, while a small subset of young 
stars are in an intermediate st ate, thought to be in "tran- 
sition" ([Skrutskie et al.lll990f ). Many of these interme- 
diate objects have optically thick "cold disks" as seen 
by large excesses at ^30 /im, but much s maller or non- 
existe nt near- and mid- infrared excesses (jBrown et al.l 
12001 . 

One model of this transition process is one where the 
inner disk is cleared of material, leaving an inner hole 
and low accretion rates, followed ~ 10 5 years later by 
a clearing of the outer disk. This inner hole could 
be caused by natu ral processes of grain growth and/or 
photo-evaporation (|Alexander et al.l l2006) , or truncation 
of th e disk due to dynamical interactions with compan- 
ions (|Lin fc Papaloizoul Il979f ). Although the most ex- 
citing interpretation for the development of gaps is ar- 
guably truncation due to planetary formation, it is also 
possible to dy namically truncate disks due to stellar com- 
panions (e.g Uensen Sz Mathieul 1X9971 iBeust fc Dutrevl 
120051: I White fc Hillenbrand 120051 ). An example of this 
is CS Cha, which wa s announced as a transitional disk 
(jEspaillat et al . 2007) s hortly after its d iscovery as a cir- 
cum binary disk llGuenther et al.1 120071 ). Some authors 
(e.g. iFurlan et al.ll2007t ) have considered the term "tran- 
sitional" to be an observational term that means the 
spectral signature of a cleared inner disk rather than an 
evolutionary term that implicitly applies to single stars. 
We use the term "cold disk" for this purpose. 

CoKu Tau/4 is a weak-lined Ml. 5 T Tauri star (El- 



Electronic address: mireland@gps.caltech.edu 

1 Data presented herein were obtained at the W.M. Keck Ob- 
servatory, which is operated as a scientific partnership among the 
California Institute of Technology, the University of California and 
the National Aeronautics and Space Administration. The Obser- 
vatory was made possible by the generous financial support of the 
W.M. Keck Foundation. 



a equivalent width 1.8 A, ICohen fc Kuhil Il979h in the 
Taurus star-forming region (l~2Myr, ~145pc). CoKu 
Tau/4 was discovered to have a large ~20-30/xm ex- 
cess and no excess at wavelengths < 8/xm through 
Spitze r Space Telescope Infrared Spectrograph observa- 
tions (|Forrest et al.ll2004l ) . The disk has been modeled as 
havin g an inner hole of radius ~10 AU fo'Alcssi o et al.l 
|2005|) , with suggestions th at this hole is due to a giant 
planet (jQuillen et al.ll2004fh 

In this paper, we describe near- infrared aperture- 
masking interferometry and imaging observations that 
demonstrate that CoKu Tau/4 is a near-equal mass bi- 
nary star system. We show that the predicted inner hole 
size from dynamical models is comparable to, but larger 
than, that predicted from radiative transfer models. Fi- 
nally we discuss whether other so-called "transitional" 
disks could be circumbinary disks. We conclude that for 
candidates in Taurus, much of the mass ratio-separation 
space for stellar companions can by ruled out by existing 
observations, but that definitively ruling out the possibil- 
ity of binarity is generally difficult for individual so-called 
"transitional" disks. 

2. OBSERVATIONS 

CoKu Tau/4 was observed with the NIRC2 camera be- 
hind Adaptive Optics (AO) at the Keck II telescope on 
2007 Nov 23 as part of an ongoing aperture-masking sur- 
vey of nearby young star-for ming associations. A perture- 
masking interferometry (e.g. lTuthill et al1l2000h is a well 
established technique for achieving the full diffraction 
limit of a single telescope, recently applied to observa- 
tions behind adaptive optics systems (e.g. iLlovd et al.l 
l2006t iKraus et al.l l2008). A 9- hole mask was placed in a 
filter wheel near a pupil plane in the NIRC2 camera, 
enabling interference fringes on 36 baselines to be si- 
multaneously recorded on the camera's imaging array. 
The observations of CoKu Tau/4 consisted of two im- 
age sets taken through a K' filter, each with eight 20 
second exposures, calibrated by two interleaved image 
sets of CX Tau. The airmass of CoKu Tau/4 observa- 
tions varied from 1.28 to 1.38, while the airmass of CX 
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Fig. 1. — Calibrated squared- visibilities for CoKu Tau 4, plotted 
against baseline projected along a position angle of 306 degrees. 
A binary model with contrast 1.25 and separation 54 mas is over- 
plotted. The dashed line shows the visibility expected for a single 
star. 
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Fig. 2. — Measured closure phases for CoKu Tau 4 as a function 
of modeled closure phases, assuming that CoKu Tau 4 is a 1.25:1 
binary with separation 54 mas at a position angle of 306 degrees. 
The dashed line shows the zero closure-phase expected for a single 
star. 

Tau (6 degrees from CoKu Tau/4) varied from 1.32 to 
1.45 over the 45 minutes it took to complete the ob- 
servations. From the Fourier transforms of the target 
and calibrator's images, squared-visibilities and closure- 
phases were extracted. The squared visibilities of CoKu 
Tau/4 were calibrated by dividing by the squared vis- 
ibilities of the calibrator, and the closure-phases were 
calibrated by subtracting the closure-phases of the cali- 
brator (which are non-zero due to instrumental effects). 
These observations were made in poor seeing conditions 
(uncorrected seeing l-2"in K-band), so our observations 
only achieved typical sensitivities of ~50:1 at the diffrac- 
tion limit, less than t he ^200:1 we achieved in Upper 
Sco (|Kraus et al.ll2008l ). These high sensitivities are pri- 
marily due to the non-redundant nature of the aperture- 
mask, which means that calibration is independent of 
wavefront structure on scales larger than a single sub- 
aperture (~lm). 

We found that only a binary model provided a good fit 
to the aperture-masking data, with a near-equal binary 
well within our detection limits. The calibrated squared 
visibilities and closure phases are shown in Figures [1] and 
[21 demonstrating that a binary solution provides an ex- 
cellent fit to the data. The binary solution based on the 
aperture-masking data is given in Table [1] 



Fig. 3. — Top line: Contour plots of the original AO image, 
only the A component (best- fit B solution subtracted), and only 
the B component (best-fit A solution subtracted) in the Br7 filter, 
demonstrating the PSF fitting technique. Bottom line: AO images 
in the K CO nt, H CO nt and Jcont filters, clearly showing the binary 
despite the low Strehl. The contours are drawn at 10% to 90% in 
intervals of 10%. 



TABLE 1 
Direct Imaging Observations 



Filter 


Wavelength 


P 


PA 


A 




Range ()im) 


(mas) 


(deg) 


(mag) 


Aperture 
K> 


Masking 








1.948-2.299 


54.1 ±0.3 


306.7 ± 0.4 


0.23 ± 0.01 


PSF Fitting 








B C ont 


1.569-1.592 


53.6 ±0.7 


306.3 ± 0.3 


0.24 ± 0.04 




2.152-2.185 


53.3 ±0.2 


307.7 ±0.7 


0.25 ±0.01 


Kcont 


2.256-2.285 


53.6 ±0.2 


307.1 ±0.2 


0.21 ± 0.01 


Bispcctrum Analysis 








J cont 


1.203-1.223 


51.7 ±0.9 


305.8 ± 0.6 


0.23 ± 0.06 


Hcont 


1.569-1.592 


53.3 ±0.2 


306.9 ±0.2 


0.21 ±0.02 




2.152-2.185 


53.9 ±0.2 


306.4 ±0.2 


0.19 ± 0.01 


Kcont 


2.256-2.285 


53.6 ±0.2 


306.0 ±0.2 


0.20 ± 0.01 



Preliminary analysis of the masking data indicated 
that the binary companion's separation and flux ratio 
would allow it to be resolved with direct imaging, so we 
obtained regular AO images on the following night under 
slightly better observing conditions (uncorrected seeing 
~l"in K-band) at an airmass of 1.06. These observations 
were obtained with NIRC2 in the narrow camera mode 
(~10 mas pix -1 ), and since the target is very bright in 
the near- infrared, we used the narrow-band filters J CO nt, 
Hcont, K cont , and Brj. Total exposure times were 40s 
in all filters except for J CO nt, where twice as many im- 
ages were taken due to the very low Strehl, for a total 
exposure time of 80 s. We show the images in Figure [3) 
Despite the low Strehl ratio and the elongation of the 
speckle halo due to a dominant wind direction, the bi- 
nary is clearly resolved in all images. 

We extracted photometry and astrometry from the di- 
rect AO imaging of these sources using the PSF recon- 
structio n technique described in iKraus fc Hillenbrand! 
(|2007bft . which iteratively fits a template PSF to the 
primary and then subtracts the secondary to fit an im- 
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proved estimate of the primary. This routine was im- 
plemented using the ALLSTAR routine in the IRAF 
package DAOPHOT (Stetson 1987), and based on previ- 
ous experience with NIRC2 imaging, we chose to model 
the PSF using a Gaussian core with Lorentzian wings. 
This method is illustrated in Figure [3] In addition to 
this method, we use d a recursive-phase algorithm (e.g 
lLohmann et al.lll983l ) to reconstruct Fourier Phase from 
bispectrum phase, and then fitted a binary model to 
the reconstructed phase of the image. In Table [TJ we 
summarize the relative astrometry and photometry that 
we measured from our direct imaging observations. The 
consistency in separation and position angle at the 1% 
level in H and K bands and the 5% level in J band de- 
spite very different PSFs and diffraction effects demon- 
strates the reliability of our results. However, the small 
reported statistical errors in Table [1] are not entirely con- 
sistent internally, requiring additional systematic errors 
to be added. Therefore, we assign a weighted mean so- 
lution of separation 53.6 ± 0.5 mas and a position angle 
of 306.4 ± 0.6 degrees. The contrast ratio in J, H and K 
is consistent with there being no significant difference in 
reddening between the components. 

3. SYSTEM PROPERTIES 

3.1. Stellar and Binary Properties 

Both components have near infrared colors that are 
consistent with other young Taurus members of spec- 
tral type ~M1-M2, suggesting that they have similar 
temperatures and that neither has a significant K band 
excess due to a hot inner disk. In order to model the 
stars and the disk, it is important to know the bolomet- 
ric luminosity of the system, and hence the reddening. 
To estimate the reddening, we used the foll o wing data: 
intrinsi c colo rs reported in Besscll fc Brettl (|1988t ) and 
iBesselll (|1990l) interpolating between giants (log(<?) ~ 1) 
and dwarfs (log(<?) ~ 5) at an assumed log(g) of 3.7, JHK 
photom etry from 2MASS, V-b and photometry averaged 
between [Cohen fe Kuhil (|1979f ) and lHanson et all ([20041 ) 
with an assumed 0.1 mag overal l uncertainty, a nd an 
Ml. 5 spectral classification from ICohen fc Kuhil (|1979l ) 
with an assumed 1 sub-class uncertainty. We arrive at 
E V -K = 2.3 ± 0.3 and Ej- K = 0.58 ± 0.04. The ra- 
tio Ev-k/Ej-k — 4.0 ± 0.6 is significantly l ower than 
the 5.6 typical of Taurus (jWhittet et al.ll200ll ). suggest- 
ing a reddening law more typical of a dense core. Based 
on Ey-K, the reddening Ay is 2.5 ± 0.3 using the rela- 
tionship A v w I.IEv^k from lWhittet et all (|2001l ). but 
given t he uncertainties in the reddening law, Ay = 3 as 
used bv lD'Alessio et all (|2005f ) is certainly consistent. 

Absolute stellar properties like masses can be diffi- 
cult to estimate because pre-main-sequence evolution- 
ary models are not well -calibrated, but based on th e 
prescription described in iKraus fc Hillenbrand! (|2007af ). 
both components are likely to have masses of ~0.5-0.6 
Mq. Using a 145 pc distance to Taurus, the system has 
an absolute K magnitude of 2.55, which implies absolute 
K magnitudes of 3.2 an d 3.4 for the c ompo nents. Ac- 
cording to the models of iBaraffe et al.l (|1998f ). assuming 
an eff ective temperat ure of ^3600 K applicable to Ml. 5 
stars (jLuhmanl Il999h . the stars are 4MYr old and are 
of mass ^0.52 and ~0.62M©. If the effective temper- 
ature is closer to 3500 K, the stars would be of slightly 
lower mass and < 3 Myr old, consistent with the canoni- 



cal age of Taurus. Th e observational scatter in the Tau- 
rus HR diagram (e.g iLuhmanl [2004) is consistent with 
CoKu Tau/4 being either single as assumed by previous 
authors or an equal brightness double. The mass ratio, 
however, should be more robust than the masses or age 
(as it is only weakly dependent on reddening) , so we can 
assert with greater confidence that Mb /Ma =0.85±0.05. 

The minimum possible semi-major axis (a) for the or- 
bit is the apparent separation divided by (1+e), which 
is 7.8/(l+e) AU based on an assumed distance of 145 pc 
for Taurus, giving an orbital period >10 years. Based 
on considerations of dynamical disk truncation in Sec- 
tion 13.21 it is likely that the system was observed near 
maximum elongation. 

3.2. The Circumbinary Disk 

We expect the inner radius of the CoKu Tau/4 disk 
to be set by tidal truncation. To simplify our discus- 
sion, we will assume that the disk is roughly co-planar 
with the binary and (as described in Section l3~Tj) that the 
bin ary components have similar masses. In this geome- 
try, lArtvmowicz fe Lubowl (|1994h calculate that the 3:1 
orbital resonance at Rd ~ 2.08a is open ed for all but im- 
proba bly small eccentricities (e > 0.03). iBeust fe Dutrevl 
(|2005f) predict a gap size of > 2.6a for eccentricities 
> 0.1 for GG Tau, which has a similar mass ratio to 
CoKu Tau/4 (and find this not to be large enough for 
that system). For very large eccentricities (e > 0.6), 
larger resonances ope n and predicted gap sizes from 
lArtvmowicz fc Lubowl (|1994h are in excess of Rd ~ 3a 
but are viscosity-dependent. The closer near-equa l spec- 
troscopic binary HD 98800B (| Akeson et alj I2007D is an 
example of a system with high eccentricity that should 
have opened large resonances in the disk: but the ra- 
diative transfer modeling has too many uncertainties to 
clearly say which resonance has opened. 

Based on the minimum semi-major axis of 
7.8/(l+e) AU from Section [2j we arrive at mini- 
mum disk truncation radius of ^16 AU if CoKu Tau/4 is 
a low (e ~ 0.1) eccentricity binary, or ~13AU if CoKu 
Tau/4 is a high (e > 0.6) eccentricity binary. These 
numbers are discrepant from the 9-12 AU truncation 
radius predicted by t he detailed spectral modeling of 
ID'Alessio et al.l ([20051 ). who examined radiative transfer 
models consisting of a geometrically thin wall and a 
range of dust types. We will examine several ways to 
resolve this discrepancy in turn: the differing radiation 
pattern of a binary versus a single star, errors in the 
stellar luminosity, and the possibility of a substantially 
different radiative transfer model. 

A binary star has a different radiation pattern to a 
single star. In particular, on a circular ring of radius Rd 
centered on the binary center of mass, flux is greatest 
at points on this ring aligned with the axis of the bi- 
nary. For an 8 AU binary (each component 4 AU from the 
center of mass) and a 14 AU ring, this geometry means 
that the flux from the binary is between 0.92 and 1.28 
times the flux from a single star with the same total lu- 
minosity. We will ass ume that the ~145 K modeled by 
ID'Alessio et al.l ([20051 ) represents the maximum temper- 
ature seen in the ring and this temperature determines 
Ls/R D - Here L$ is the total stellar luminosity and Rd 
is the dust radius from the center of mass. Therefore, a 
28% increase in flux F could represent a ~13% increase 
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in the ring radius from radiative transfer modeling. The 
temperature contrast around the ring for optically thin 
dust with opacity k\ oc A -1 goes as F 1 / 5 , meaning that 
the temperature around the ring ranges from ~136K to 
145 K. Gray dust would vary from -133 K to 145 K. Al- 
though this contrast is only moderate, it may present an 
additional challenge in detailed spectral fitting, as the 
models of D' Alessio et al.l (|2005|) required emission from 
a single temperature only. 

One way to reconcile the gap radius from radiative 
transfer modeling with that from dynamical considera- 
tions would be if the luminosity of CoKu Tau/4 is sig- 
nificantly underestimated. If we keep the dust temper- 
ature and therefore Ls/R 2 D constant, increasing L$ in- 
creases the model Rd- An uncertainty in dist ance due 
to the depth of Taurus (e.g. iTorres et'afll2007f ) can not 
reconcile this discrepancy, as both Rd determined from 
radiative transfer modeling and the minimum Rd deter- 
mined by dynamical considerations both increase linearly 
with the distance estimate. However, an uncertainty in 
the reddening of CoKu Tau/4 translates to a luminosity 
uncertainty. The probable redde ning calculated in Sec - 
tion l3.1l is lower than that used bv lD'Alessio et al. ( 2005), 
but uncertainties are considerable, and it may be possi- 
ble that a larger bolometric reddening could increase the 
gap radius derived from radiative transfer modeling. 

In order to independently estimate the dust radius 
from radiative transfer modeling and appreciate possi- 
ble further complexities, we will examine a very simple 
model. This model contains small grains of amorphous 
olivine (i.e. the limit where the grain radius is much 
smaller than the blackbody emission pea k), with opti- 
cal constants from iDorschner et alj (|1995h . We assign a 
dust tempera ture Tp=145K a nd a stellar luminosity of 
0.61 -L© as in iD'Alessio et al\ (j2005| ). an d a stellar tem - 
perature of an Ml. 5 star of T s =3600 K (|Luhmanlll999l ). 
The stellar radius for each of the two stars is 3.7 Rq, 
equivalent to a single star of radius 5.2 Rq. We assume a 
blackbody spectrum for the star 2 . The dust is assumed 
optically thin (which could represent the outer layers of 
an optically-thick disk), so that the dust radius is then 
given by: 



Rd = Rs\ 



1 jK X B x (T s )d\ 
2 / Kx B x (T D )dX 



(1) 



In this simple model, the dust radius is Rd =15 AU, 
roughly consistent with disk truncation radii from dy- 
namical considerations. With the many possible parame- 
ters to tweak in a more complete model, and the possibil- 
ity of the optically-thin dust being additionally heated by 
surrounding dust, increasing Tr> and driving the model 
radius outwards, it appears likely that a detailed radia- 
tive transfer model could indeed be made to match the 
dynamical truncation radius. The edge in the spectrum 
of disk emission from CoKu Tau/4 at — 9 /im coincides 
with a sharp rise in silicate emissivity, meaning that the 
interplay between temperature, extinction, geometry and 
dust composition is all the more complex for this source. 
We hope that this preliminary discussion will motivate 



other authors to examine possible models in more detail. 



4. DISCUSSION AND CONCLUSIONS 

The circumbinary nature of the CoKu Tau/4 disk begs 
the question: are other so-called transition disks also 
likely circumbinary in nature? The other well-known 
cold disks in Taurus, DM Tau, LkCa 15, UX Tau and 
GM Aur, are also generally assumed to be single stars. 
Although all these systems differ from one another, CoKu 
Tau/4 is perhaps most unique because of its very low 
accretion rate. Our aperture-masking observations of 
these systems (paper in preparation) can so far elim- 
inate stellar (mass ratio q > 0.1) companions over a 
20-160 mas projected separation range, but can say lit- 
tle about the possibility of closer companions. Accurate 
(er < lkms' 1 ) radial velocity monitoring of these sys- 
tems is required to determine if any of them harbor a 
close (<4AU) companion. Radial velocity monitoring 
is most applicable to DM Tau, where the cleared inner 
hole is expected to have on ly a 3 AU radius (|Calvet et al.l 
120051 ). iSimon et alj (|2000[ ) used resolved CO 2-1 obser- 
vations of disks around several T Tauri stars to measure 
their masses directly. These measurements would have 
resulted in the total mass of the stellar primaries and any 
unresolved companions. The total masses of 0.84 ± 0.05, 
0.97 ± 0.03 and 0.55 ± 0.03 for GM Aur, LkCa 15 and 
DM Tau respectively provide convincing evidence that 
these systems do not harbor companions with mass ra- 
tios q > 0.3. 

In summary, CoKu Tau/4 is a binary star system with 
a current projected separation of ~7.8AU, and a mass 
ratio near unity. The disk surrounding CoKu Tau/4 is 
therefore a circumbinary disk, with the disk inner edge 
set by dynamical truncation. The dynamical disk trun- 
cation radius of > 13 — 16 AU is inconsistent with the 
disk radius derived from previous spectral modeling, but 
it is likely that a different choice of reddening law and 
correction, dust composition and geometry will resolve 
this small discrepancy. We suggest that systematic radial 
velocity monitoring and aperture masking surveys are re- 
quired to determine if other so-called transition disks are 
in fact multiple systems. 
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2 Assuming instea d a NextGen model spectrum from 
IHauschildt et ail [fL999h at 3600 K and log(g)=4.0 reduces R D by 
2AU, but has far too much TiO absorption when compared to an 
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